(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organization 

International Bureau 




llllllllllllllllllllllllllllllllllllllllllllllllllli 



(43) International Publication Date (10) International Publication Number 

30 January 2003 (30.01.2003) pCT WO 03/007809 A2 



(51) International Patent Classification^: A61B 5/00 

(21) International Application Number: PCT/CA02/01081 

(22) International Filing Date: 16 July 2002 (16.07.2002) 
(25) Filing Language: English 



(81) Designated States (national): AE, AG, AL, AM, AT, AU, 

AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CO, CR, CU, 

CZ, DE, DK, DM, DZ, EC, EE, ES, H, GB, GD, GE, GH, 
GM, HR, HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, 
LK, LR, LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, 
MX, MZ, NO, NZ, OM, PH, PL, PT, RO, RU, SD, SE, SG, 
SI, SK, SL, TJ, TM, TN, TR, TT, TZ, UA, UG, US, UZ, 
VN, YU, ZA, ZM, ZW. 



(26) Publication Language: 



(30) Priority Data: 

60/305,092 
09/985,436 



16 July 2001 (16.07.2001) 
2 November 2001 (02.1 1.2001) 



(71) Applicant (for all designated States except US): ART, 
ADVANCED RESEARCH TECHNOLOGIES INC. 

[CA/CA]; 2300 Alfred-Nobel Blvd., Saint-Laurent, 
Quebec H4S 2A4 (CA). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): HALL, David, 
Jonathan [GB/CA]; 3445 me Drummond, app. 1005, 
Montreal, Quebec H3G 1X9 (CA). KHOL-BAREIS, 
Mathias [DE/DE]; Drosselveg 21, 53489 Sinzig (DE). 



English (84) Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZM, ZW), 
Eurasian patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), 
European patent (AT, BE, BG, CH, CY, CZ, DE, DK, EE, 
ES, FI, FR, GB, GR, IE, IT, LU, MC, NL, PT, SE, SK, 
TR), OAPI patent (BF, BJ, CF, CG, CI, CM, GA, GN, GQ, 
GW, ML, MR, NE, SN, TD, TG). 



US 
US 



Published: 

— without international search report and to be republished 
upon receipt of that report 

For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



(74) Agents: OGILVY Renault et al.; 1981 McGill College 
Avenue, Suite 1600, Montreal, Quebec H3A 2Y3 (CA). 



(54) Title: CHOICE OF WAVELENGTHS FOR MULTTWAVELENGTH OPTICAL IMAGING 



00 



O 



0.012 
0.01 



^ 0.008 



< ^ 



0.006 



B 0.004 
0.002 
0 



oxy-Hb 10 pM 



deoxy-Hb 10 pM 

water 

fat 




650 700 750 600 850 900 
wavelength [nm] 



950 



(57) Abstract: The present invention 
relates to a method for wavelength 
selection in a multi-wavelength 
TPSF-based optical imaging system. 
This consists of identifying several 
chromophores in a highly turbid 
medium and selecting optimized 
wavelengths whereby using these 
wavelengths optimizes the deduction of 
the chromophore concentrations. Such 
chromophore concentrations may be 
combined to deduce other properties of 
the turbid medium. 
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CHOICE OF WAVELENGTHS FOR MULTIWAVELENGTH 

OPTICAL IMAGING 

Field of the Invention 

The present invention relates to the field of optical imagin'g In which objects 
5 which diffuse light, such as some human body tissues, are imaged using signals 
resulting from the injection of light into the object and detection of the diffusion of 
the light in the object at a number of positions. More particularly, the present 
invention relates to the choice of wavelengths for multiwavelength optical 
imaging in order to provide enhanced information. 

10 Background of the invention 

Time-domain optical medical images show great promise as a technique for 
imaging breast tissue, as well as the brain and other body parts. The objective is 
to analyze the temporal point spread function (TPSF) of an injected pulse as it is 
diffused in the tissue, and the information extracted from the TPSF is used in 
15 constructing a medically useful image. 

For example, one can extract time-gated attenuation information from the TPSF 
which improves the image spatial resolution over previous continuous wave 
methods. However, it is unclear whether such improvements in image spatial 
resolution are adequate for diagnosing breast cancer based on morphology. 

20 An alternative approach is to use the TPSF to decouple the light attenuation into 
absorption and scattering components. This extra information, which cannot be 
obtained from continuous wave methods, may be clinically useful. Moreover, 
one can obtain the tissue absorption spectrum by performing time-domain 
measurements at multiple wavelengths. In tissue there are several molecules 

25 which absorb the light and are known as chromophores. Spectroscopic analysis 
of the tissue absorption spectrum permits chromophore concentrations to be 
measured. Furthermore, combination of the chromophore concentrations can 
yield physiological information, as opposed to morphologic information, which 
could provide a more medically useful image. 

30 The problem is one of knowing which are the dominant chromophores to include 
in a tissue model and then choosing the "best" wavelengths to deduce their 
concentrations most accurately. 
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Summary of the invention 

It is an object of the invention to improve image quality in TPSF-based optical 
images by choosing an efficient combination of wavelengths and combining 
information from the combination of wavelengths. 

5 It is an object of the present invention to provide an objective method for 
choosing the wavelengths for a multiwavelength TPSF-based optical imaging 
approach. For a given set of chromophores, the best selection of the 
wavelengths is performed for the set as a whole as opposed to choosing the 
best wavelength for each chromophore individually. Furthermore, hardware 
10 constraints can be taken into consideration in order to optimize the selection of 
wavelengths for a given device. 

Brief description of the drawings 

Fig. 1 illustrates the absorption spectra used of oxy-Hb, deoxy-Hb, pure water 
and lipid; 

15 Fig. 2 illustrates the inverse of the condition number of the hemoglobin specific 
absorption matrix as a function of wavelength X^ being plotted for a system of a) 
two wavelengths where the other wavelength is fixed at X2 = 850 nm, b) three 
wavelengths where the other wavelengths are fixed X2 = 850 nm and ;\.3 = 758 
nm and c) four wavelengths where the other wavelengths are fixed X2 = 850 nm, 

20 X3 = 758 nm and 7A = 800 nm; 

Fig. 3 illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb as a function of XI and X2. The plot is 
symmetric with respect to the diagonal. Regions of high values indicate 
combinations of wavelengths advantageous for spectroscopy; 

25 Fig. 4 illustrates the deviation of calculated saturation and true saturation 
(S(calc)-S(true)) for a model tissue containing 15 ^xM [HbT], S(true) = 25%, 50% 
and 75% and a lipid concentration of 40%. Two wavelengths at 760 and 850 nm 
were used to fit [oxy-Hb] and [deoxy-Hb]. The sensitivity with respect to wrong 
assumptions of lipid and water concentrations are shown; 

30 Fig. 5A illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb and lipid for a fixed wavelength X = 830 nm as 
a function of A.1 and X2, The islands of high values indicate advantageous 
wavelengths (scaling 0-0.01); 
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Fig. 5B illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb and lipid for a fixed wavelength X3 = 830 nm as 
a function of XI and X2 (same as Fig. 5A but scaling 0-0,0005); 

Fig. 6A illustrates the inverse of condition number C for the specific absorption 
5 spectra of oxy-Hb and deoxy-Hb, lipid and water for two fixed wavelength at X3 = 
760 nm and X4 = 830 nm as function of and X2 (scaling 0-0.0015). Regions 
of high values are advantageous for spectroscopy; 

Fig. 6B illustrates the inverse of condition number C for the specific absorption 
spectra of oxy-Hb and deoxy-Hb, lipid and water for two fixed wavelength at X3 = 
10 760 nm and X4 = 830 nm as function of X^ and X2. (same as Fig. 6A but scaling 
0-0.0005); 

Fig. 7 illustrates the estimation of deviations from true saturation values for a 
model tissue of [HbT] = 20 |j,M, S= 75%, a lipid concentration of 40% and true 
water concentration corresponding to 0-100% water. Three wavelengths at 760, 
15 780 and 850 nm were used for back calculation of S, shown here as a function 
of assumed water concentration; 

Fig. 8 illustrates the estimation of the influence of errors (noise) in |Lia on the 
calculated Hb concentrations and saturation values. A model p-a-spectrum based 
on 20 |laM [HbT], S=50%, and a lipid and water concentration of 30% and 40% 
20 was assumed. Matrix inversion was performed for wavelengths 760, 790, 830 
and 850 nm. plotted is the change in calculated [oxy-Hb], [deoxy-Hb] and 
saturation value when the ij,a value at a single wavelength was changed by 
+0.0001 mm"\ This plot suggests that noise at 830 nm translates in the highest 
noise in saturation values; 

25 Fig. 9 illustrates the estimation of the recovery of saturation values based on 
different wavelength combinations. A model tissue of 20 juM [HbT], a true 
saturation of S=75%, lipid and water concentration of 40% were used, in the 
lower plot an offset of 0.0005 mm"^ independent of wavelength was added to the 
model tissue ^la-spectrum (no offset in the upper plot). Plotted are deviations of 

30 the saturation values due to matrix inversion and the true 75% value. The 
following wavelength combinations were used: 1) 760 nm and 850 nm, 2) 760, 
830 and 850 nm, 3) 760, 780, 830 and 850 nm, 4) 750-850nm, 5) 720-850 nm, 
6) 720-900nm; 

Fig. 10 illustrates the estimation of the recovery of saturation values based on 
35 different wavelength combinations. A model tissue of 20 p,M [HbT], a true 
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saturation of S=50%, lipid and water concentration of 40% were used, in the 
lower plot an offset of 0.0005 mm"^ independent of wavelength was added to the 
model tissue jua-spectrum (no offset in the upper plot). Plotted are deviations of 
the saturation values due to matrix inversion and the true 75% value. The 
5 following wavelength combinations were used: 1) 760 nm and 850 nm, 2) 760, 
830 and 850 nm, 3) 760, 780, 830 and 850 nm, 4) 750-850nm, 5) 720-850 nm, 
6) 720-900nm; 

Detailed Description of the invention 

In accordance with the present invention, there is provided a method for 
10 selecting wavelengths for multlwavelength optical imaging. 

Tissue Chromophores 

The dominant near infrared chromophores contained in breast tissue are 
considered to be hemoglobin (Hb) in its oxygenated (oxy-Hb) and deoxygenated 
(deoxy-Hb) forms, water and lipids. Fig. 1 shows the absorption spectra of oxy- 

15 Hb (at lOiLxM concentration), deoxy-Hb (at lOjxM concentration), pure water 
(100% concentration), lipid (absorption spectrum of olive oil has been used to 
estimate the absorption spectrum of fat). There are other interesting near 
infrared chromophores, such as glucose and cytochrome c oxidase, but their 
absorption contribution in the breast is considered negligible compared to the 

20 aforementioned chromophores. 

Physiological information 

Potentially useful physiological information about the breast tissue can be 
obtained from concentrations, [], of the chromophores. The total hemoglobin 
concentration, [HbT], defined as [HbT] = [oxy-Hb] + [deoxy-Hb], is related to the 

25 local vascular density. Since cancer is commonly associated with an increase in 
vascularisation (angiogenesis), a measurement of [HbT] could be medically 
useful. The fraction of hemoglobin that binds to oxygen is known as the oxygen 
saturation, S, and defined as S=[oxy-Hb]/[HbT]. Increased metabolic activity 
increases oxygen demands which decreases the oxygen saturation. Since 

30 cancer is commonly associated with increased metabolic activity, a 
measurement of S could also be medically useful. 

Wavelength Choice 

Historically as the biomedical optics field evolved the wavelengths were chosen 
for each chromophore individually by observing strong near infrared spectral 
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features for the given chromophore and using the closest hardware-available 
wavelength. Many researchers also used the isobestic wavelength of oxy-Hb 
and deoxy-Hb, the wavelength where their absorption per concentration are 
equal, since this wavelength is insensitive to the oxygenation state of the 
5 hemoglobin and can be related to the [HbT]. 

However, the question both posed and addressed here is that for a given set of 
chronnophores what are the optimal wavelengths to use in order to deduce the 
concentration of each chromophore? It is interesting to note that the isobestic 
wavelength used by many researchers turns out not to be one of the 
1 0 wavelengths of choice. 

It is an object of the present invention to provide an objective method for 
choosing the wavelengths for a multiwavelength TPSF-based optical imaging 
approach. For a given set of chromophores, the best selection of the 
wavelengths is performed for the set as a whole as opposed to choosing the 

15 best wavelength for each chromophore individually. Moreover, it is also possible 
to investigate scenarios such as the influence on determining chromophore 
concentrations under certain assumptions about the concentration(s) of other 
chromophore(s) in the set. Furthermore, hardware constraints can also be taken 
into consideration in order to optimize the selection of wavelengths for a given 

20 device. Fortunately, the recent advent of turn-key, pulsed, tunable near infrared 
wavelength lasers has permitted more viable availability of near infrared 
wavelengths. 

Experimental brute force approach 

One possible approach to optimize the choice of wavelengths for a given set of 
25 chromophores is to conduct a brute force experimental study. This would consist 
of performing numerous experiments where different combinations of 
wavelengths are evaluated for the given set of chromophores at known 
concentrations until the optimum combination for deducing their concentrations 
is found. Obviously, this approach is likely to be highly time-consuming and it is 
30 not always trivial to provide a set of chromophores at known concentrations, 
particularly in the case of in vivo breast tissue. 

Matrix inversion sensitivity approacli 

An alternative approach which avoids the numerous experiments of the 
experimental brute force approach is a matrix inversion sensitivity approach. 
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The equation which needs to be solved can be written for each wavelength as: 

i 



5 , Ma(.^) = J^rn„M3)-^i 

i 

where is the measured absorption coefficient, nia is the specific absorption 
coefficient of the different chronnophores and 0\ is the corresponding 
concentration. 

10 This is written is nnatrix fornn as: 

Ha = Wl-C 

where printing in bold indicates a matrix or vector. |Lia is a vector with a number of 
rows corresponding to the number of wavelengths (n^,). c is a vector with the 
number of rows corresponding to the number of chromophores (nc). M is a 
1 5 rectangular matrix of size n^ x nc- 

If n;^ = nc the system can be solved by matrix inversion c=M"^jia and if n^ > nc the 
system is overdetermined and can be solved by the pseudo-inverse M* = (M^ 
M)"'' where is the transposed matrix of M. 

c = (M'*)Kla 

20 The pseudoinverse M* is an n^ x nc array which is unique. If M is square (i.e. not 
overdetermined) the = For given (i.e. chosen) wavelengths the 

pseudoinverse M'*' can be precalculated once and the matrix inversion 
corresponds to a simple matrix multiplication. This is the basis for the calculation 
of chromophore concentration. 

25 One means to quantify the expected sensitivity of a matrix inversion of a matrix 
M with respect to small errors in the data is the condition number C which is 
defined as: 

C = norm (M) • norm (M"^) 

C gives an indication of the accuracy of the results and is an estimate of the 
30 cross-talk between the different channels (i.e. chromophores concentrations). 
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Values of C near 1 indicate a well-conditioned matrix, large values indicate an ill- 
conditioned matrix. The condition number is closely related to singular value 
decomposition (SVD) as it is the ratio of the largest and the smallest singular 
value of a matrix. 

5 The matrix M for oxy-Hb and deoxy-HB at A. = 760 and 770 nm is 

^0.3871 0A465A = 760nm 
■^0.3280 0.1625^ = 770/27^ 

A matrix inversion is possible as the rank (M) = 2, however the absorption at the 
two wavelengths is 'similar'. The condition number is C = 20.49. Choosing the 
wavelengths to be A. = 760 and 850 nm gives the matrix 

0.3871 0.1465 

10 M = 

0.1729 0.2645 

Inspection by eye already shows that the absorption is very 'different'. This is 
confirmed by the condition number: C=3.206. In what follows below the inverse 
of the condition number is plotted and analyzed. It has value between 0 and 1. 
1/C close to 1 means 'orthogonal' spectra and low sensitivity to cross-talk. Small 
15 values of 1/C mean an ill-conditioned matrix. To find the best wavelengths, 1/C 
is calculated as a function of a wavelength. The wavelengths that give the 
highest values of 1/C are the best for a calculation of chromophore 
concentrations and the subsequent physiological information such as oxygen 
saturation, S. 

20 Model absorption spectra were generated with the absorption spectra of Fig. 1 
based on estimations of [HbT], S, lipid and water concentration. Matrix inversion 
based on different sets of wavelengths were performed to recover these 
parameters. These parameters were compared with the true ones for the 
different wavelengths and the sensitivity to noise or measurement offsets 

25 considered. 

Assuming that we fit for the hemoglobin concentrations only and assuming 
certain values for water and lipid concentration, for a x-wavelengths matrix 
inversion, the best combination of wavelengths to give a well-conditioned matrix, 
the sensitivity of calculated values of oxy-Hb and deoxy-Hb concentration and 
30 oxygen saturation for variations of lipid or water concentration and sensitivity of 
S to measurement noise have been determined. 
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In Fig. 2 the inverse of the condition number is shown for matrices of oxy-Hb and 
deoxy-Hb specific absorption coefficients for 2, 3 and 4 wavelengths. In each 
case one wavelength (A.i) was varied between 650 and 950 nm while the 
remaining wavelengths were fixed X2 = 850 nm (2-wavelength system), = 850 
5 nm and A.3 = 758 nm (3-wavelength system), and X2 = 850 nm, X3 = 758 nm and 
%A = 800 nm (4-waveiength system). Fig. 2 indicates that the selection of two 
wavelength at X\ = 850 nm and A,2 = 700 nm gives the highest values of 1/C and 
when the wavelength range is restricted via hardware constraints to > 750 nm, a 
system that includes the peak wavelength of deoxy-Hb close to 760 nm is 
10 advantageous. It does not matter whether two or more wavelengths are used. 
This somewhat counterintuitive result is valid only, without measurement noise 
and noise in the background absorption. 

Fig. 3 further highlights this finding for a two-wavelength matrix inversion. In this 
figure 1/C is plotted as a function of both at X\ and Xz in the range 650-950 nm. 
15 The plot is symmetric with respect to the diagonal. Regions of high 1/C-values 
can be chosen and the corresponding 'good' wavelengths can be read off the 
axis, it is apparent that (with the restriction to > 750 nm) the one wavelength 
should be close to 760 nm while the other one can be within the range 830-900 
nm without substantially affecting the condition number. 

20 Using the spectra shown in Fig.1, model tissue absorption spectra were 
generated. Based on matrix inversion values of [oxy-Hb], [deoxy-Hb] and S were 
backcalculated and the sensitivity to incorrect assumptions about the [water] and 
[lipid] tested. One approach is to take the measured |Lia spectra and subtract 
water and lipid absorption corresponding to certain assumed concentrations. For 

25 the data shown in Fig. 4, a model tissue containing 15 inM [HbT], (true) 
saturation values of S = 25%, 50% and 75% was used. Lipid concentration was 
40%. It was tested how a misjudgement of water concentration affects the 
recalculated S value. To test the error in a simple two-wavelengths-fit (760 and 
850 nm), the assumed lipid concentration was varied between 0 and 100%. 

30 When the assumed water concentration is right (lower three lines in Fig. 4), the 
deviation in saturation between true and calculated values is < +2% (obviously 
with zero error for the right lipid concentration of 40%). A misjudgement about 
the water concentration by 20% (upper lines in Fig. 5) results in additional errors 
in S(calc)-S(true) of up to 2% for S=75%, 4% for S=50% and 8% for S=25%. 

35 These errors in S are a function of the underlying tissue absorption coefficients. 
The values here give an indication about the order of magnitude. 
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Having a system with more than two fit-parameters, best wavelength 
combinations, for a three-components system of oxy-Hb, deoxy-Hb and lipid 
system, for a four-components system of oxy-Hb, deoxy-Hb, lipid and water, and 
the sensitivity of calculation of S to noise at the different wavelengths have been 
5 determined. 

In Figs. 5A and 5B the inverse of the condition number is plotted for a three 
wavelengths system based on the oxy-Hb, deoxy-Hb and lipid specific 
absorption spectra as a function of Xi and A,2. The third wavelength was fixed at 
X3 = 830 nm. Again, the plot is symmetric with respect to the diagonal. From Fig. 
10 5A it is apparent, that there are three "islands" of high 1/C values. Unfortunately, 
all of these island would include wavelengths outside an imposed hardware 
constrained wavelength range of 750 to 850 nm. Plotting the same data in a 
different scale (Fig. 5B) shows that there is just a single preferential combination 
within this hardware constrained wavelength range: 760 and 780 nm. 

15 Equivalent to Figs. 5A and 5B, the inverse of C for a 4-wavelengths system is 
plotted in Fig. 6A and 6B. Again, the difference between them is the scaling. 
Two wavelengths were fixed at A.3 = 760 nm and X4 = 830 nm. Including the 
wavelengths outside the 750-850 nm range there appear four preferential 
combinations. Restricting the wavelength range to 750-850 nm there are just two 

20 advantageous region (marked by the white rectangle In Fig. 6B): 780 nm and 
850 nm, and 780 nm and 815 nm. 

From the analysis based on matrix condition numbers, the best wavelength 
combinations for 2, 3 and 4 wavelengths measurements are the following: 



wo 03/007809 



PCT/CA02/01081 



-10- 



Table 1 

Wavelength combinations for 2, 3 and 4 wavelengths measurements 





Wavelength 


Best wavelengths (nm) 




fit for 


see 




Range 


x^ 


X2 


X3 


X4 




Fig. 




650-950 nm 


700 


>860 






Oxy-Hb, deoxy-Hb 


2,3 




750-850 nm 


760 


850 










3-X 


650-950 nm 


700-760 


830 


925 




+ lipid 


5A,B 






830 


860-870 


925 










750-850 nm 


760 


780 


830 










650-950 nm 


760 


830 


860 


925 


+ lipid, + water 


6A,B 






700 


760 


830 


925 








750-850 nm 


760 


780 


830 


850 


Best combination 








760 


780 


815 


830 







5 Furthermore, it must be pointed out that including more wavelengths does not 
increase the condition number E.g. for the four chromophores and all 
wavelengths in the range 750-850 nm, 1/C = 0.000314. This is lower than the 
value (1/C = 0.00036, compare with Fig. 6B) when only four wavelengths (760, 
780, 830 and 850 nm) are used. In a system without noise and no other 
10 chromophores than the four considered here, a 4-wavelengths system is the 
optimal. 

While certainly only a 4-wavelengths measurement allows [oxy-Hb], [deoxy-Hb], 
[lipid] and [water] to be determined, and a 2-wavelengths system (see Fig. 5) is 
not sufficient, the question is posed whether a 3-wavelengths measurement 

15 might supply S values with a high enough precisipn. In this case the 
concentration of one chromophore (water or lipid) must be guessed and the 
corresponding absorption subtracted from the measured jua values. This was 
tested with a model absorption spectrum and is shown in Fig. 7 for 760, 780 and 
850 nm. True water concentration was varied between 10 and 100% (the 

20 different lines), and the difference between calculated and true saturation values 
plotted as a function of assumed water concentration. For instance, for a true 
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water concentration of 50%, a misjudgnnent of the water concentration by 10% 
results in an error in S by about 4%. 

Up to now, only "perfect" data sets were considered with no noise. In real 
situations there are problems due to measurement noise that is random for the 
5 different wavelengths; unknown chromophores in the tissue, i.e. there is a 
background absorption coefficient the spectrum of which we do not know; and 
possible systematic errors in the primary |LXa recovery. 

There are an ample number of parameters which can be considered and as 
examples the following two questions are considered. First, is the oxygen 
10 saturation more susceptible to noise at certain wavelengths? Second, what is 
the influence of an offset in the |Lia data? 

The influence of errors (noise) in |Lia on the calculated Hb concentrations and 
saturation values was estimated in a model tissue based on 20|liM [HbT], S= 
50% and a lipid and water concentration of 30% and 40% respectively. Matrix 

15 inversion was performed on the lUa - spectrum of this model tissue for 
wavelengths 760, 780, 830 and 850 nm. The sensitivity to noise (i.e. variations 
in i^a) at the different wavelengths was estimated by changing the absorption 
coefficient at a single wavelength by +0.0001 mrn\ In Fig. 8, the change in 
calculated [oxy-Hb], [deoxy-Hb] and oxygen saturation value due to this "noise" 

20 is plotted. This figure shows that the change in oxygen saturation value is about 
-2% for changes at 760 nm, <0.5% at 780nm, while it translates to a variation of 
+6% at 830 nm. 

There are two conditions that produce an offset in the measured |ia-spectra with 
respect to the true values. First, the algorithm for jua-calculation based on TPSF- 

25 based optical imaging might lead to a systematic offset e.g. due to residual 
crosstalk between absorption and scattering parameters. Second, the tissue 
absorption might have a background of unknown origin (chromophore). Under 
both conditions the fitting of ^a-data with the four chromophores is hampered. 
The effect of such an offset for different wavelength combinations is estimated 

30 with a model spectrum of 20 ^iM [HbT], S=75% and water and lipid concentration 
of 40%. An offset of 0.0005 mm"^ was added to the ^la-values independent of 
wavelength. The effect on the calculated oxygen saturation values is known in 
Fig. 9 for combination of 2, 3 and 4 wavelengths as well as continuous spectra 
between 750-850, 720-850 and 720-900 nm. It is apparent that the lowest error 

35 in S is achieved by the 4-wavelengths combination. Including more wavelengths 
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increases the error. In Fig. 10 the same calculation was done, however, for a 
true oxygen saturation value of 50%. Here the lowest error is achieved by the 
720-850 nm wavelength range, while using less wavelengths or increasing the 
fitting range to 900 nm results in larger errors. 

5 Based on the assumption that the dominant tissue chromophores are oxy-Hb, 
deoxy-Hb, water and lipid and analysis of the matrix condition number, 
measurements at the wavelengths 760, 780, 830 and 850 nm supply an optimal 
data set when the wavelength range is limited to 750-850 nm under ideal 
conditions. As shown in Figs. 5 and 6, inclusion of shorter and longer 

10 wavelengths promise a better matrix inversion. Under real conditions there is no 
clear-cut answer about the improvement when more wavelengths are included 
(see Figs. 9 and 10). It might be advantageous to reduce measurement noise at 
4 wavelengths due to longer scan times rather than to include more 
wavelengths. As demonstrated in Fig. 8, to achieve an optimal accuracy the 

15 noise level at different wavelengths has to be adjusted which might require 
different measurement times at certain wavelengths. 

Strictly speaking the work presented here was achieved by optimizing a 2- 
wavelength system and then optimizing a 4-wavelength system where 2 of the 
wavelengths were fixed at the optimized 2-wavelength system values. Whilst this 
20 is easier to display graphically, preferably all 4 wavelengths would be permitted 
to vary in a global optimization process. Fortunately, for the specific example 
presented here when all 4 wavelengths are permitted to vary the same optimal 
solution is found. However, this may not be true for all situations and a global 
optimization is preferred. 

25 It is also understood that it will be obvious to those skilled in the art that the 
same approach for choosing optimal wavelengths can be applied to optical 
absorption spectroscopy in general. For example, in other embodiments of the 
present invention the method of the present invention is also used for choosing 
the optimal wavelengths for analyzing the components of paints, pharmaceutical 

30 products, food, grain or any other turbid media. 

It is also understood that the proposed method applies both to the analysis of 
absolute chromophore concentrations as to their changes or relative 
concentrations. 

While the invention has been described in connection with specific embodiments 
35 thereof, it will be understood that it is capable of further modifications and this 
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application is intended to cover any variations, uses, or adaptations of tlie 
invention following, in general, the principles of the invention and including such 
departures from the present disclosures as come within known or customary 
practice within the art to which the invention pertains and as may be applied to 
5 the essential features herein before set forth, and^as follows in the scope of the 
appended claims. 
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CLAIMS 

1. A method of optical imaging of turbid media using a plurality of 
discrete wavelengths in an optical imaging system, the method comprising the 
steps of: 

selecting a set of chromophores for characterizing a property of the 
turbid media; 

defining parameters of the system including at least a number of said 
discrete wavelengths, a value of each of said wavelengths, source power and 
detector aperture for each of said wavelengths, a choice of image algorithm and 
source/detector geometries, a choice of source and detector and noise 
characteristics; 

fixing a value of all of said parameters except a plurality of said 
parameters values to be optimized; 

determining an optimal value for each of said parameter values to be 
optimized as a function of a performance of the system in measuring a 
concentration of said chromophores in said turbid media for characterizing said 
property as a whole; and 

using said optimal value for each of said parameter values in imaging 
said turbid media. 

2. The method of claim 1, wherein said imaging is medical imaging, said 
highly turbid medium being body tissue and said property is physiological. 

3. The method of claim 1 or 2, wherein said parameter values to be 
optimized comprise a value of each of said wavelengths. 

4. The method of claim 3, wherein said parameter values to be optimized 
further comprise said number of said discrete wavelengths. 

5. The method of claim 4, wherein said step of determining comprises 
fixing said number of discrete wavelengths at each of a plurality of numbers, and 
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determining an optimized performance of the system in measuring a 
concentration of said chromophores in said turbid media at each of said plurality 
of wavelengths, and selecting one of said plurality of numbers having a best 
optimized performance. 

6. The method of claim 3, wherein said step of determining an optimal 
value for each of said parameters comprising minimizing a condition number of 
a matrix of specific absorption coefficients of said chromophores as a function of 
wavelength, 

7. The method of claim 4, wherein said step of determining an optimal 
value for each of said parameters comprises minimizing a condition number of a 
matrix of specific absorption coefficients of said chromophores as a function of 
wavelength. 

8. The method of claim 5, wherein said step of determining an optimal 
value for each of said parameters comprises minimizing a condition number of a 
matrix of specific absorption coefficients of said chromophores as a function of 
wavelength, 

9. The method of claim 1, wherein said step of determining comprises 
empirically determining said performance of the system for a range of said 
values for each of said parameter values to be optimized. 

10. The method of claim 2, wherein said plurality of chromophores 
comprise at least oxy-hemoglobin and deoxy-hemoglobin. 

11. The method of claim 10, wherein said chromophores are water, lipids, 
oxy-hemoglobin and deoxy-hemoglobin. 



12. The method of claim 10, wherein said body tissue is breast tissue. 
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IS. The method of claim 10, wherein said number of wavelengths 
selected is from 2 to 4. 

14. The method of claim 13, wherein said number is 4. 

15. The method of claim 11, wherein values of said wavelengths are 760 
nm, 780 nm, 830 nm and 850 nm. 

16. The method of claim 1, wherein the step of determining an optimal 
value of said parameters to be optimized comprises: 

deriving an inherent wavelength-dependent sensitivity to noise in 
calculating said chromophore concentrations, and 

determining an optimal correlation of said sensitivity and at least one 
other of said parameters, 

17. The method of claim 16, wherein one of said parameters to be 
optimized is a distribution of an acquisition time at each of said wavelengths. 

18. The method of claim 1, wherein one of said parameters to be 
optimized is a distribution of an acquisition time at each of said wavelengths. 

19. The method of claim 18, further comprising a step of determining a 
minimum value for said acquisition time at which said performance of said 
system attains a minimum threshold value. 

20. The method of claim 1, wherein one of said parameters to be 
optimized is at least one of said source power and said detector aperture for 
each of said wavelengths. 
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21. The method of claim 20, further comprising a step of determining a 

minimum value for an acquisition time at which said performance of said system 
attains a minimum threshold value. 
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